Detecting Josephson effect in the excitonic condensate by coherent emission of light 
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Coherent emission of light by a split excitonic Bose-Einstein condensate - excitonic Josephson 
junction - can be extremely sensitive to the Josephson phase established across the junction. As 
a result of this, the emission can be redistributed between different directions and even cancelled. 
The cases of the dipole- and of the quadrupole-active excitons are considered. 
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Bose-Einstein condensation of low density excitons attract much of experimental and theoretical efforts jjj-ji). 
Depending on the nature of the excitons, the excitonic condensate may display a variety of properties specific for 
superfluids and superconductors. One of these is Josephson effect. In Ref. [|| it has been shown that an optical analog 
— i, of the Josephson effect - optical Josephson effect - can be observed in the excitonic condensate. This effect can be 
revealed by the phase-sensitive transmition of light || . In Ref. it was suggested that the Josephson current can 
be detected in the case of the condensate of indirect excitons due to spatial separation of holes and electrons. 

Reports on realization of the excitonic condensate |2|,[| are based on observing a phonon assisted excitonic photo- 
luminescence whose spectrum is proportional to the excitonic distribution in the energy space. Such a method does 
not test directly coherence of the excitonic cloud. In Ref. || the excitonic photoluminescence from the quantum 
wells subjected to magnetic field have exhibited large noise which was considered as an evidence of the excitonic 
condensation stimulated by the field and characterized by finite size domains. Recently, several proposals have been 
discussed (8 11 . These rely on observing light produced by the condensate of optically active excitons. In Ref. ||, the 



two photon emission has been shown to be a unique property of the condensate. Statistical properties of the emitted 
radiation have been addressed in Refs. P,|lol- The polarization of the coherent emission from the two component 
excitonic condensate has been discussed very recently [jllj . It has been shown that it is sensitive to the relative phase 
q ■ between the components. 

In this paper it is pointed out that the excitonic Josephson effect between two confined excitonic condensates may 
affect significantly the intensity of the light emitted by such condensates. Thus, observing this effect can be a direct 
' proof of the excitonic condensation. Besides that, various properties of the excitonic Josephson junction can be tested 
, by means of measuring the emission intensity. 

As long as excitons can transform into light, the local polarization is proportional to the excitonic operator \E'. In 
t-H | the excitonic condensate the mean tp = becomes finite due to a formation of the off-diagonal long range order 
Accordingly, the mean electric polarization (P) ~ (W) = ijj exp(— iojt) becomes finite, where u> is the frequency of 
the excitonic transition. This implies that the emission becomes coherent, as opposed to the case when only thermal 
excitons emit f[fj[| . 

Here a possibility of observing a dependence of the emission intensity on the relative phase ip between two excitonic 
condensates will be analyzed. A simple general analysis shows that this dependence should exist. Indeed, let us 
assume that two identical confined condensates are spatially separated by some displacement Xo and do not overlap 
significantly. Then, the two-condensate wavefunction can be represented as 



V(x) =Vi(x)+e iV ^(x-x ), (1) 

C 

where V'i( x ) an d e l¥> ?/>i( x — x o) are the wave functions of the first and the second condensates, respectively. The 
electric and magnetic fields of the coherent radiation at the excitonic frequency ui in the far zone are given by the 
polarization density which is proportional to ip. This leads to the differential intensity W(q) ~ IV'ql 2 °f emission of 
light with the wave vector q, where ipq stands for the spatial Fourier component of V>( x )- Thus, from the representation 
rS | (Q) one finds 



W(q)=2Wi(q)(l + cos(^-qx )), (2) 

where Wi(q) denotes the intensity due to one condensate only when the second is not present at all. Eq.(||) shows 
that the differential emission intensity of two condensates depends on their relative phase tp. 

If xq is much smaller than the wavelength A = 2ir/q of the emitted light, one may neglect the term qx in Eq.(||). 
Then, the total intensity W^ tot ^ = J dnW(q) as an integral over the 47r-solid angle (n = q/q) will be very sensitive 
to the relative phase. Specifically, 

W {tot) = 2Wf°' ) (l + cos^) (3) 
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where w[ tot) = fdnWi(q). It can be seen that for tp — ±n the emission can be practically cancelled. This is an 
example of the dark state which can be realized in the excitonic Josephson junction. 

However, if the spatial separation xq of the condensates is comparable with or larger than A, the inegration over 
the total solid angle may suppress the dependence on the phase tp in Eq.fQ). Nevertheless, as it will be seen below, 
the total intensity of emission by two condensates of optically active excitons can remain strongly sensitive to the 
phase for large xq. In contrast, in the case of the quadrupole-active excitons in C112O, the dependence on the phase 
practically disappears for xq >> A. 

First, let us analyze the case of optically active (OA) excitons confined inside, e.g., two identical spatially separated 
cubes with the side 2a. The excitonic condensate is characterized by the density n = \ipi\ 2 which is taken the same 
for both condensates. It is naturally to consider a limit when the excitonic healing length I = l/y/na s is much 
smaller than the wavelength A of the emitted radiation. Indeed, for typical achievable densities n w 10 18 cm~ 3 [QJJ] 
and the excitonic scattering length a s taken as a typical excitonic radius i=a 10~ 7 cm, one finds I ~ 10~ 6 cm. Under 
this condition the condensate density can be taken as a constant inside the cubes. To specify a relative position of 
the cubes, these are assumed to be identically oriented and displaced by x > 2a along one of their 4th-fold symmetry 
axis, which will be called the x-axis. The excitonic polarization vector P is directed along another 4th-fold symmetry 
axis, which will be called the z-axis. Thus, the polarization density vector (0, 0, P z ) can be represented as jf3| 

P z = d$(0)*, (4) 

where d denotes the interband dipole matrix element; ^ stands for the operator of the excitonic center of mass 
motion, and $(r) is the wave function of the relative motion of the electron and the hole forming the exciton. In the 
condensate phase at the temperature T = the operator 'J is to be replaced by the condensate wave function ip [[l4| . 
Given (||) for a pure condensate case, one finds the differential intensity Wi(q) in Eq.(||) as 

w l(q) = wo ( M°*)™fa)Mw) Y (5) 

\ q sm 9 cos 9 sin tp cos q> J 

where within the orthogonal coordinate system x,y,z the angles 9 and 4> give the direction of the vector q so that 
the components of this vector are q z — qcos9, q x = q sin 6 cos <f> and q y = q sin 9 sin <fr; the introduced constant is 
Wq — 8c|d<I>(0)| 2 n/7r, where c stands for the speed of light. A simple analysis of the angular dependence (||) shows 
that, when qa » 1, the emission pattern is characterized by a very anisotropic angular distribution. Almost all the 
radiation occurs inside 4 narrows cones of the angular area S « l/(qa) 2 << 1. These cones are collinear with the x and 
y axes. Thus, employing Eq.(||) in Eq.(|J), one can introduce four integrated intensities Wx tot \ W^"^, Wy tot \ W^°^ 
emitted into the corresponding directions x > 0, x < 0, y > 0, y < 0. It is important to note that the intensities 
Wy tot ^ — emitted perpendicular to the vector x , are practically insensitive to the phase ip. In contrast, the 

intensities Wx tot \ W^° l \ emitted collinear with x , are very sensitive to the phase tp. Specifically, one finds 

w (tot) = 2W -M)(i + CQS ^ _ Xoq ^ ; = 2W^ tot) {l + cos{p + log)), (6) 

w (tot) = w (tot) = 327rc | d $(o)| 2 n a 2 sin 2 (oa) (7) 

in the limit qa >> 1. Note that for x q = pn where p is integer, the emission in the both directions x > 0, x < can 
be completely suppressed by choosing tp = or tp = ir depending on parity of p. For xoq =/= pit, the emission can be 
cancelled in either x > or x < direction. It is important to note that such a phase dependence of the emission is 
a direct consequence of the excitonic condensation. No such an effect will occur when no condensate is present. 

The phase difference can be induced by a difference of the chemical potential created by, e.g., some external field 
or due to a slight difference in the exciton number. The condensates can also be prepared initially with a certain 
phase difference imposed by a resonant light ||. If the two condensates are connected by a constriction, an excitonic 
flux from one condensate to another will pass through this narrow bridge and exhibit a sudden change of the velocity. 
The velocity of the condensate is associated with the gradient of the phase tp. Therefore, the velocity change in the 
constriction will result in a jump of tp between two condensates. As far as the total emission intensity (g) is very 
sensitive to the phase, the phase dynamics can be studied by means of measuring the total emission intensity. 

A special consideration should be given to the case of excitons in C1L2O. The paraexcitons have been reported 
to form the condensate 0. In Ref. jj] the condensation of orthoexcitons has been reported as well. However, the 
orthoexcitons are characterized by fast Auger decay Thus, a possibility of achieving a coherence time long 

enough to observe the Josephson effect in a short living condensate of the orthoexcitons is questionable. In what 
follows the case of the long living paraexcitons will be considered. In the non-stressed crystal no direct radiation can 
be observed from the paraexcitons. The phonon assisted emission is inherently incoherent. However, if a mechanical 
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stress producing the deformation Uij is applied, a quadrupole and magnetic dipole radiation can be detected from the 
radiatively recombining paraexcitons ||l^,|l7|. Below it will be shown that, if the geometry of the confinement of the 
two paraexcitonic condnesates is properly chosen, the excitonic Josephson effect can be revealed in this case as well. 
For sake of simplicity, the cubic crystal of CU2O will be approximated by isotropic medium. Then, instead of Eq.(^), 
one may write for the electric polarization density 

Pi = 5^Qii$(o)a,-* (8) 

where Qij is the interband quadrupole tensor of the optical transition. In isotropic medium the second-rank tensor Qij 
can be induced by the deformation Uij so that Qij — guij in the limit of small deformations, where g stands for some 
constant, which can be determined from the experimental dependence of the intensity of the direct recombination 
emission as a function of the applied stress fl7|| . Let the deformation be along the z axis, so that the only non- vanishing 
component is u zz . Thus, again, as in the case of the OA excitons, the polarization is along z-axis (P z — gu zz <fr(Q)d z il)). 
Now, however, this polarization is concentrated very close to the boundaries where the gradients of the excitonic 
condensate are the largest. Because of this, the emitted radiation becomes more isotropic than in the case of the 
OA excitons considered above. Consequently, the term qx under the cos in Eq.(|J) will practically erase the phase 
dependence unless xq is less than q = 2tt/\ (so that Eq.(^) can be employed). Thus, the considered above geometry 
of the two cubes is not appropriate for the paraexcitons. Instead, it is reasonable to consider a different arrangement. 
Specifically, let us consider a slab whose width along the z-axis is 4zo << A, and whose base is a square with the 
side 2a >> A oriented along the x, y axes. This slab is divided in half by the plane perpendicular to the z-axis. The 
Josephson phase ip is assumed to exist between these two halves. Given the representation (||), calculations of the 
emission intensity in this case show that the factor W\ in Eq.(||) is to be replaced by 

w (p) ( q ) = w (p) ( sm(aq x )sm(aq y )sm(z q z ) \ 2 
1 ° \ sin sin cos / 

where Wo = (16/7t)c|Q zz $(0)| 2 ti and again q z — qcos9, q x — qsm.6coS(f> and q y = g sin sin </!>. It can be seen that 
even in the limit qa » 1 the radiation is no more concentrated in the plane 9 = tt/2. Instead, it is spread with 
respect to the axial angle 9 around the directions 6 = tt/4, 37r/4. However, the radiation is still very anisotropic 
with respect to the polar angle (p. It is equally emitted along the directions cf> = 0, tt, ±7r/2 into a very small angle 
S(j) ~ « 1- Let us calculate the integrated intensity emitted into one of these directions (e.g. = 0) 

by means of integrating Eq.(||) over Jj^ dcos9 J**^ d(f> , with the replacement W^ p \q) — > Wi(q) performed. This 
yields 

W {p) = 16Trcnq 3 z%a\Q zz <f>(0)\ 2 (l + cos<^) (10) 

in the limit zq << A. 

Now let us discuss a role of the thermal component. The intensities (||), (0), ([!(]) of the coherent emission are not 
proportional to the total volume of the condensates. In the case of the OA excitons the intensity (||), (R) scales as 
~ a 2 and for the paraexcitons the intensity ( |Tc| ) is ~ a in the limit a — > 00. In contrast to this, the intensity of the 
incoherent emission due to the thermal excitons should be proportional to the total volume occupied by the 

excitons. Specifically, W^ m ' ~ a 3 for the OA excitons, and W^ m ^ ~ zqo 2 in the case of the paraexcitons. Thus, one 
may expect a masking of the coherent emission by the thermal component. 

However, it should be noted that, while the incoherent radiation is almost isotropic, the coherent emission is very 
concentrated in certain directions. Therefore, a comparison should be performed between the intensities integrated 
over the small solid angles, to where the most of the coherent radiation is emitted. In the case of the OA condensate, 
the integrated incoherent emission intensity is ~ a? /(qa) 2 ~ a, where it was taken into account that the coherent 
emission is concentrated into the cones of the angular area 1=3 l/(qa) 2 . Accordingly, the ratio of the incoherent to 
the coherent integrated intensities is ~ 1/a — > in the limit a — > 00. Thus, the coherent radiation should always be 
resolved for large enough values of a. 

A different situation occurs for the paraexcitons. The angular area of the coherent emission "cones" is 5(f) ~ 1/ qa. 
Thus, the integrated incoherent intensity is ~ a 2 j qa ~ a, and the ratio of it to the intensity (|l^) is idependent of the 
largest size a of the system. In what follows the intensity of the incoherent emission will be calculated for the case of 
the paraexcitons, and the condition insuring that the coherent effect can be resolved will be derived. 

In order to find the incoherent intensity, one should employ the representation (^|) for the polarization. At T 7^ 0, 
the operator $ = tp + tp' contains the condensate part ip and the normal part ip' . The differential intensity of the 
incoherent emission at the wave vector q is W tyln \(\} ~ (V'qV-'q) where the averaging is performed over the thermal 
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ensemble. The average can be calculated by employing the Bogolubov method H]. Specifically, the thermal part is 
represented as 

ip' = ^2{u k a k +v* k al), (11) 
k 

where u k and v k stand for the Bogolubov amplitudes: a k destroys and at creates a normal excitation at the state k 
characterized by the energy E k . Here the quasi-particles can be characterized by the quasi- momentum q inside the 
slab. Thus, the index k must be replaced by q in Eq. (|il"|) ■ A direct conversion of the excitation with given q into a 
photon with almost the same q (in the limit qa >> 1) is considered. In the contribution due to the normal component, 
the only non-zero means are (et q a q ) which are the population factors of the quasi-particles. Thus, calculating the 
emission intensity into a unit solid angle due to the thermal paraexciton, one finds 

^ (m) (q) = ^IQ,^(0)l^sin 2 (2g) g6(l 7' 2 + ^ |2) , (12) 
62-k exp{E q /l ) — 1 

where V = IQzqo 2 stands for the total volume of the excitonic cloud; the quantum depiction term ~ |v g | 2 , which is 
finite at T = 0, has been neglected. In Eq.(|l^) the amplitudes and the energy are H] 

\u q \ 2 = £q+9 ; x ^ + E \ \v q ?= eq + 9 Z£- E \ E q = ^e 2 q + 2g ex ne q . (13) 

In here e q = q 2 /2m denotes kinetic energy of a free exciton, and g ex — 4ira s /m stands for the two-body excitonic 
interaction constant expressed in terms of the scattering length a s (in units % = 1 ). In Eq. ([l2]) the main term 
proportional to the total volume of the system is selected. In the limit n\ 3 >> 1 and for the gas parameter na s << 1, 
one can safely replace the Bose factor in Eq.([l2|) by T/E q unless T is extremely low. Here it is assumed that T 
is of the order of the temperature T c » 3.3n 2 / 3 /m of the excitonic Bose- Einstein condensation (for one component 
excitons) for typical densities ~ 10 18 cm~ 3 . Under these circumstances the expression ([l2| ) is weakly sensitive to the 
interaction constant g ex in the limit E q <<T. Specifically, it can be seen that W^ n > changes by only a factor of 2, 
while g ex varies from zero to such values that the excitonic spectrum becomes sound like E q — y/ g ex n/m q. In the 
latter situation the g-dependent factor in Eq.([l2"|) becomes q 6 T/2e q , that is independent of g ex . In the opposite limit 
(g e x — ► 0), this factor is q 6 T/e q . 

The total intensity of the incoherent radiation emitted into the narrow solid angle S(f>, to where the most part of 
the coherent radiation is collected, is » W^ ln \(\} / qa. Thus, the ratio of the integrated incoherent intensity to that 
given by Eq.([To|) (for tp = 0) is 

The coherent effect can be well resolved if £ < 1. It can be seen that for not very small zq, one obtains £ << 1 as 
long as T < T c . Indeed, making an estimate of £ for the typical values of the wavelength A « 600 nm, the density 
n = 10 18 cm~ 3 , the temperature T = T c /2 and 2itzq/\ = 1/3, one finds £ w 10~ 2 . Let us also estimate the absolute 
intensity of the coherent emission from the paraexcitonic condensate. To this end, the following maximum value of Q zz 
is chosen Q zz w ea 2 B , where e stands for the electron charge and as — 0.5A denotes Bohr's radius. ^(O)! 2 = l/ira 3 x 
p3| is expressed in terms of the excitonic radius a ex w 10A. Thus, for the chosen above values, one arrives at the 
intensity ( ffO ) as ~ 10 _5 W for a/A = 10 3 . It is worth noting that the intensity of the coherent emission from 

the condensates of the OA excitons is a factor of ~ d 2 a/q 3 z 2 Q 2 z larger than the above estimate. Depending on the 
choice of d which can be varied by external fields jl8j , this factor can enhance the above estimate by several orders of 
magnitude. 

In summary, the coherent radiation from two excitonic condensates can be sensitive to their relative phase. In the 
case of the OA excitons, the emission can be totally suppressed or significantly redistributed in the opposite directions 
by means of imposing appropriate values of the phase. In the case of the paraexcitons in CuiO under mechanical 
stress, the dependence of the emission on the phase can be also very distinct, if the two condensates form a double 
layer whose thickness is less than the wavelength of the emitted radiation, and if the mechanical stress is applied 
perpendicular to the plane of the layer. 
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